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Abstract. We investigate perturbations of a class of spherically symmetric solutions in mas¬ 
sive gravity and bi-gravity. The background equations of motion for the particular class of 
solutions we are interested in reduce to a set of the Einstein equations with a cosmological 
constant. Thus, the solutions in this class include all the spherically symmetric solutions 
in general relativity, such as the Friedmann-Lemaitre-Robertson-Walker solution and the 
Schwarzschild (-de Sitter) solution, though the one-parameter family of two parameters of 
the theory admits such a class of solutions. We find that the equations of motion for the 
perturbations of this class of solutions also reduce to the perturbed Einstein equations at 
first and second order. Therefore, the stability of the solutions coincides with that of the 
corresponding solutions in general relativity. In particular, these solutions do not suffer from 
non-linear instabilities which often appear in the other cosmological solutions in massive 
gravity and bi-gravity. 


Contents 


1 Introduction 1 

2 Review of bi-gravity 3 

3 Bi-spherically symmetric background solutions 5 

4 Linear perturbations 7 

5 Second-order perturbations 8 

6 Conclusions and discussion 9 

A Concrete examples of background solutions 10 

A.l Bi-cosmological solutions 10 

A.2 Bi-Schwarzschild de Sitter solutions 11 

B Bianchi identity 12 

C Additional gauge symmetry of linear perturbations 13 

D Hamiltonian analysis of linear perturbations 15 

D.l Odd mode perturbations 15 

D.2 Even mode perturbations 18 


1 Introduction 

Massive gravity is one of the potent candidates for modified theory of general relativity. As 
early as in 1939, a linear theory of massive gravity was proposed by Fierz and Pauli (FP) [1]. 
In order to avoid the inconsistency on massless limit of this theory [2, 3], a nonlinear extension 
of the FP theory was considered [4]. Boulware and Deser, however, found that the nonlinear 
theory simply extended from the FP theory contains an unphysical ghost degree of freedom 
(BD ghost) [5]. Because of this ghost problem, a healthy theory of non-linear massive gravity 
had not been established for a long time. 

Recently, de Rham, Gabadadze, and Tolley (dRGT) proposed a mass potential which 
can remove the BD ghost mode in a decoupling limit [ 6 , 7], and Hassan and Rosen have 
hnally proven that the dRGT massive gravity theory is free from the BD ghost without 
taking the decoupling limit [8-10]. The dRGT theory of massive gravity has three parameters, 
graviton mass m and coupling constants of nonlinear self interactions, 013 , 0 : 4 . Complementary 
approaches of the BD ghost problem are studied in refs. [11-14]. 

Massive gravity includes a non-dynamical tensor field called fiducial metric, /^i^, in 
order to construct a mass potential. For example, the FP and the original dRGT theories 
are constructed by adopting the Minkowski fiducial metric. Hassan and Rosen proposed an 
extended theory of dRGT massive gravity with a hducial metric being dynamical as well by 
introducing the Einstein-Hilbert term of the hducial metric in the action. They proved that 
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this theory is also free from the BD ghost [15]. Since this theory contains two symmetric 
dynamical tensor fields of metrics, it is called bi-gravity theory. 

The tests of massive gravity and bi-gravity using cosmological and black hole solu¬ 
tions have been explored intensively. In dRGT massive gravity, several types of exact, 
homogeneous, and isotropic solutions have been known so far. One example is the open 
Friedmann-Lemaitre-Robertson-Walker (FLRW) solution with the flat fiducial metric in the 
FLRW slice [16]. The second-order perturbations of this solution show nonlinear instability 
and hence this solution is not viable unfortunately [17-19]. It should be noted that similar 
solutions with an anisotropic fiducial metric are known to be stable [20-22]. Another example 
of cosmological solutions is that with a fiducial metric which is flat but expressed in terms of 
nontrivial coordinates. This class of solutions can be divided into two types. The first type 
includes the solutions found in refs. [23-25], which exist for the whole parameter region of as 
and 04 , while the other type includes the solutions found in refs. [26, 27], which exists only 
for a one parameter family in the parameter space ( 03 , 04 ). Though the perturbations of the 
former type of solutions have already been studied in refs. [28-31], the perturbations of the 
latter type of solutions have not yet been investigated. Therefore, in this paper, we focus 
on the latter type of solutions. It should be noted that cosmological solutions with non-flat 
fiducial metrics are also studied in refs. [32-35]. 

The situations on cosmological solutions in bi-gravity are similar. The cosmological 
solutions found thus far are divided into two classes [36-39] . The first class is a solution with 
diagonal metric tensors [40-42] . The perturbations of this class of solutions have been studied 
intensively [43-55]. On the other hand, for the other class of a solution with off-diagonal 
components of physical or fiducial metric tensor, the perturbations have not yet been studied, 
similarly to the case of massive gravity. It should be noticed that coupling between matter 
and (bi-)metrics is a nontrivial issue in bi-gravity and is studied in refs. [56-63]. 

A lot of static and spherically symmetric solutions have also been found up to now. 
The classification of such spherically symmetric solutions is studied in refs. [39, 64, 65]. The 
exact Schwarzschild (-de Sitter) solutions are classified to the following three classes. The 
first class is a solution with diagonal metric tensors [66-69] and linear perturbations of this 
class of solutions are studied in ref. [69-73] in the framework of both massive gravity and 
bi-gravity. The second class of solutions is a solution with a off-diagonal metric tensor and 
arbitrary 03 and 04 [74], and the perturbation of this solution is studied in ref. [69, 75]. The 
last class is a solution with a off-diagonal metric tensor and a special choice of the parameters 
03 and 04 [76-78], where linear perturbations have been studied only in massive gravity with 
a flat fiducial metric [79] and not in bi-gravity. 

In the present article, we will give a unihed and general treatment for solutions with 
an off-diagonal metric tensor in massive gravity and bi-gravity belonging to a one parameter 
family of 03 and 04 , which include both cosmological [26, 27] and spherically symmetric 
black hole solutions [64, 76-79]. We will hnd that the equations of motion for this class 
of solutions exactly reduce to those of general relativity (GR) with a cosmological constant 
not only at the background and linear (first-order) perturbation level but also at the level of 
quadratic (second-order) perturbations. This result shows that massive gravity and bi-gravity 
can allow any spherically symmetric solution of GR including its stability, the evolution of 
linear perturbations, and the backreaction from linear perturbations, while it simultaneously 
implies that one cannot distinguish massive gravity or bi-gravity from GR by using spherically 
symmetric solutions and their perturbations at least up to quadratic order. 

Our paper is organized as follows. In the next section, we briefly review the theory 
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of bi-gravity (and massive gravity as a trivial case of a fixed fiducial metric) and derive 
the equations of motion in a general setting. In section 3, we derive a generic non-diagonal 
spherically symmetric background solution. Then, we investigate linear perturbations around 
those background solutions in section 4. There we will see that the terms coming from the 
mass potential must vanish in order to satisfy the Bianchi identity. In section 5, we investigate 
higher order perturbations and find that the same results as the linear perturbations apply 
for the quadratic perturbations. The final section is devoted to summary and discussion. 
Some details will be given in the appendices. 


2 Review of bi-gravity 

In this section, we give a brief review of bi-gravity. Bi-gravity is a theory consisting of 
two dynamical tensor fields, g^i, and called physical and fiducial metrics, respectively. 
Massive gravity can be understood as a special case where the fiducial metric is fixed and 
non-dynamical. Its action is given by the Einstein-Hilbert term for each metric with the 
interaction term Smass and matter actions: 

S = I d^xV^R[g] + I d^xy^Rif] 

T'S’massiS)/] T 'S'matteriS'] T 'S'matter [/]) (2-1) 

where i2[-] is the Ricci scalar and k represents the ratio of the effective Planck masses for g^^ 
and In the case of k = 0, the tensor field does not have its kinetic term and hence 
is non-dynamical. This case corresponds to the massive gravity theory originally proposed 
by de Rham, Gabadadze, and Tolley [6, 7]. S'matterb] and 5matter[/] are the matter actions 
coupled to g and /, respectively, 

‘S'matter [g] = j d^xy/^ 

■^matter’ (^■^) 

5^atter[/] = J (2.3) 

Here we implicitly assume the matter actions possess the two general covariance with respect 
to gfj,y and ffj,y separately, though the full theory does not have such a symmetry. It should 
be noted that another type of matter coupling, which does not possess the two general 
covariance, is also studied by [56-63]. The energy-momentum tensors coming from 5matterb] 
and Smatterl/j are dehned as 

rp fi _ ^ lip <^>5'matter [g] 

rr p _ ^ ^5'matter[/] ^ 

Due to the two general covariance of matter actions we assumed, both energy-momentum 
tensors are conserved, that is, = 0 and = 0, where and are 

the covariant derivatives with respect to g^y and respectively. Hereafter, we will omit 
the suffixes g and / when no confusion is expected. 


(2.4) 

(2.5) 
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Now the interaction term Smass is tuned to be free from the BD ghost mode and given 


by 


Af^ r 

Smassb,/] =2 <YxY~g2riY'^Piei{'y), 

i=0 

(2.6) 

where i-th order contributions 64 ( 7 ) are given by 


60 ( 7 ) = 1, 

61 ( 7 ) = Tr[ 7 ], 

(2.7) 

(2.8) 

62 ( 7 ) = \ (Ti'[ 7]2 - Tr[72]) , 

(2.9) 

63 ( 7 ) = ^ (Tr[7]^ - 3Tr[7]Tr[72] 2Tr[73]) , 

(2.10) 

64 ( 7 ) = det( 7 ). 

(2.11) 

with 


7^ = 

(2.12) 


and m, /3j being free parameters of the interaction term. Since they always appear in the 
combination mPPi, essentially there are hve free parameters. The space of parameters corre¬ 
sponds to the one of the three parameters of the dRGT theory, m, as, 04 , and two cosmological 


constants, , and their relations are given by^ 

m^/3o = — A^^^ -|- m ^(6 -|- das + 04 ), (2-13) 

m^/3i = m^(—3 — 3as — a 4 ), (2-14) 

m^/32 = m^(l-I-2as-I-a 4 ), (2-15) 

m^/3s = m^(—as — a4), (2-16) 

m^/34 = +m^ai. (2-17) 


Taking the variation of the action with respect to and we will obtain the 
equations of motion for the two tensor fields. The equations of motion for are given by 

G[9l% + (2.18) 

where G[g]^^u is the Einstein tensor constructed from g^^ and 

= 2w? , (2.19) 

+/33 (e2(7)7't - eiY){YYu + (7^)''i2)] • (2-20) 

^ It is useful to rewrite the action in term of these parameters as follows, 

4 

^ UiBiilC) + , 

■i=2 

with = 6^u — 7^iv and a 2 = 1. 


‘S'mass [^5/] — 


Ml, 


/' 
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The indices here are raised or lowered by The equations of motion for are given by 

+ ( 2 . 21 ) 

where G[f]^i, is the Einstein tensor constructed from f^i, and 

= -^sgn(det7) + . (2.22) 

The indices here are raised or lowered by 


3 Bi-spherically symmetric background solutions 


Here, we attempt to classify some of the spherically symmetric solutions in bi-gravity and 
identify those which obey the same equations of motion as in general relativity. These classes 
of solutions include the cosmological and black hole solutions known so far [26, 36-39, 64, 79]. 
Let us consider the following bi-spherically symmetric metrics: 


g^ydx^dx'^ = gQQ{t,r)dt^ + 2gQi{t,r)dtdr + gii{t,r)dr^ + R{t,r)'^d^‘^, 
f^ydx^^dx’' = foo{t,r)dt‘^ + 2foi{t,r)dtdr + fu{t,r)dr'^ + A^{t,r)R^{t,r)dn‘^, 

with dQ"^ = dO^ + sin^ OdcjP'. The matrix g~^f takes the following form. 


{-g-^fru = 


(rVJi 0 0 \ 

{r^fYo (5■V)^ 0 0 

0 0 A‘^{t,r) 0 

Vo 0 0 A‘^{t,r)J 


(3.1) 

(3.2) 


(3.3) 


and from these ansatz it is straightforward to see that the square root of the above matrix 
is of the form 





/a{t,r) b{t,r) 0 ^ \ 

c{t,r) d{t,r) 0 0 

0 0 A{t,r) 0 

\ 0 0 0 A{t,r)J 


(3.4) 


It should be emphasized that the following discussion does not rely on the concrete expressions 

of a{t,r), b{t,r), c{t,r), and d{t,r), but rather relies only on the fact that is of the form 

of eq. (3.4). 

As explained earlier, we are interested in the case where the equations of motion for 
both metrics reduce to the Einstein equations with cosmological constants at the background 
level. Therefore, in order for X(^g-^^y to be a cosmological term, the non-trivial off-diagonal 
components, 

= -m^S - 2 A + {A-3){A- 1)03 + {A - (3.5) 

X(g)^o = -m^c[3 - 2 A + (A- 3}(A - Ijas + (A - ifai], (3.6) 

must vanish. We focus on the case of b{t,r) / 0 or c{t,r) / 0, and A / 1, since with 
b{t, r) = c{t, r) = 0, we will obtain a diagonal metrics as mentioned in Sec. 1 and the 
perturbations of such diagonal solutions have already been studied. 
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For non-diagonal solutions we are interested in, the condition that eqs. (3.5) and (3.6) 
vanish leads to 


A{t,r) = 2«3 + «4 + 1± Va 3 + « 3 -a 4 + l ^ 

Q!3 -|- 0(4 

Another requirement necessary for to be a cosmological term is 

X(,)% - = (1 - A)C{t, r)[A-2 + {A- l)a 3 ] = 0, 

where, we have defined C{t,r) as 

2 — A{a + d) + ad — be 


C{t, r) = 


(1 - A )2 


(3.7) 


(3.8) 


(3.9) 


With C{t, r) = 0, at least three eigenvalues of 7 ^ 1 / are equal to A. This class of solutions 
includes the cosmological solutions found in refs. [23-25] and the Schwarzschild solutions 
obtained in ref. [75]. The perturbations of those solutions have already been studied in detail 
in refs. [28-31]. 

In the present study, we therefore concentrate on the case with C(t, r) ^ 0. In this case, 
the solution to eq. (3.8) is 


.4 = I±^, 

1 +as 


(3.10) 


Here we have assumed that as 7 ^ —1. Equations (3.7) and (3.10) are consistent provided 
that the parameters of the theory, as and 0 ( 4 , satisfy 


0 — 1 -|- 0:3 + — 04 {— 02 ~ ■ 


(3.11) 


This is equivalent to the condition that the two branches of the solution (3.7) degenerate. 
Thus, we see that only the particular one-parameter family of as and 04 satisfying (3.11) 
admits the class of solutions we are focusing on. Note that when eq. (3.11) is fulfilled A 
can also be expressed simply as A = —132/Ids- Note also that the range of 04 is limited as 
04 = (03 + 1/2)2 ^ 3/4 > 3/4^ 

In this one-parameter family of 03 and 04 with eq. (3.11), the interaction terms for bi- 
spherically symmetric metrics (3.1) and (3.2) with eq. (3.10) are of the form of a cosmological 
term. For the interaction term gives 


X, 

^( 9 ) ^ ~ ^ 


with 


= m2(A-l)+A(^^), 


while for 


( f) 

^eff ~ sgn{ad — be) 


m? A — 1 



(3.12) 


(3.13) 


A 


(3.14) 

(3.15) 







In the case of dRGT massive gravity, f^i, is not a dynamical but a fixed metric, and hence 
we need not consider the equations of motion for f^i,. 

The class of solutions with C{t,r) / 0 includes the cosmological solutions [26, 27], 
the black hole solutions [76-79], the Lemaitre-Tolman-Bondi (LTB) solution [27], and the 
Reissner-Nordstrom (RN) solution [77]. Since the equations of motion for (and, in fact, 
those for as well) reduce to the Einstein equations with a cosmological constant, any 
spherically symmetric solution in GR is also a solution of the one-parameter subclass (3.11) 
of bi-gravity and massive gravity with a suitable fiducial metric. In appendix A, we present 
some examples of bi-FLRW and bi-Schwarzschild-de Sitter solutions belonging to this class. 


4 Linear perturbations 


Now, we analyze linear perturbations around bi-spherically symmetric solutions given in the 
previous section. The two tensor fields of metrics are perturbed as 

dfiv ~ 9^iv T Sg^ui 

fliv = Uu + 

The first-order perturbation, 5'y^iy, of 7 ^ 1 / is defined as 

\/ g~^f = + O(second-order perturbations), (4.3) 

which can be written in terms of the metric perturbations by solving the following equations. 


(4.1) 

(4.2) 




(4.4) 


where dg^i, = g^^dgpi, and Sf^i, = g^^6fpu- For our purpose we do not need the explicit 
form of the solution to the above equation, though it is obtained for a general fiducial metric 
in ref. [80]. Actually, without the explicit form of (^ 7 ^^, we can directly calculate from 
eq. (4.3) with eq. (3.4) as 




C{t,r) 


/O 0 0 

0 0 0 
0 0 57^3 

\0 0 -<57^2 


0 \ 
0 

-<57^3 

(57^2 / 


(4.5) 


Since the Einstein tensor satisfies the Bianchi identity X^p^G^y[g] = 0 and the energy- 
momentum tensor is conserved, the tensor Xf^g'^^y also satisfies Xp,X^g'^^y = 0. As demon¬ 
strated in appendix B, this requirement leads to a stronger condition 


= 0, (4.6) 

which yields = 0 for a, 6 = 2 ,3. (Note that we are interested in the case with C{t, r) / 0.) 
Since 

eq. (4.6) also implies 5X(^j^^y = 0. Thus, the equations of motion for the linear perturbations 
5gpy and 5fpy reduce to the linearized Einstein equations. 
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In order to see the implications of the equation (4.6) in more detail, we express 
in terms of the metric perturbations. Since only the angular components of enter the 
equation (4.6), we only have to deal with the angular components of the equation (4.4), which 
can easily be solved because 7 “;, = A5°'b for a,b = 2, 3. In fact, eq. (4.4) reduces to 

2Ad^\ = -A^dg\ + Sf\ = 0. (4.8) 

To sum up, the equations of motion for the first-order perturbations are equivalent to 
the following three equations; 



(4.9) 


(4.10) 

A^5gab - 5fab = 0 . 

(4.11) 


This is one of the main results of this investigation. The equations of motion for the pertur¬ 
bations of the two metrics coincide with the perturbed Einstein equations, though and 
^f^iu are subject to eq. (4.11). 

Then let us count the number of graviton degrees of freedom for this perturbed system. 
Each symmetric tensor field of metric has ten components, and there are, respectively, four 
constraints (the Hamiltonian and momentum constraints) in eqs. (4.9) and (4.10), since those 
equations are the same as the perturbed Einstein equations. Eurthermore, eq. (4.11) gives 
three constraints among the angular components of the perturbed metrics. We have four 
spacetime coordinates and hence there are four gauge degrees of freedom representing the 
choice of coordinates. In addition to those familiar gauge degrees of freedom, it turns out 
that there still remains another gauge transformation retaining the equations of motion (4.9), 

(4.10) , and (4.11), as explicitly shown in appendix C. Note that this gauge degree of free¬ 
dom corresponds to the ambiguity of the linear perturbations mentioned in ref. [79] for the 
Schwarzschild-de Sitter solution in the dRGT theory. Thus, the number of the remaining 
degrees of freedom is 10x2 —4x2 — 3 — (4-t-l) = 4, which coincides with that of two massless 
gravitons. We can confirm that this is consistent with the result of the Hamiltonian analysis 
given in appendix D: there are ten first class constraints and twelve second class constraints, 
and hence there are 8 (= 40 — 10 x 2 — 12) degrees of freedom in phase space. 

The above analysis can be applied to dRGT massive gravity only with eqs. (4.9) and 

(4.11) because the derivations of these equations do not depend on the equation of motion 
for Since, in this case, is composed of Stiickelberg fields, the condition (4.11) just 
determines perturbations of Stiickelberg fields. The remaining variables Sg^^ are governed by 
the Einstein equations and additional gauge symmetry appears as gauge degree of freedom 
for Stiickelberg fields. 

5 Second-order perturbations 

In the previous section, we have shown that the first-order perturbations obey the perturbed 
Einstein equations and hence the behavior of the perturbations coincides with that of GR, 
though 5g^y and are subject to eq. (4.11). One may then ask the question as to how 
one can discriminate this class of solutions in bi-gravity from the corresponding solutions in 
GR. 
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One possibility is to take into account the back reaction on the physical metric g^i, from 
at second order. For this purpose, we incorporate second-order perturbations as follows: 


9^v — ^9fiv 9^ fj-vi (^- 1 ) 

= Z/ii/ + Sffiv + (5-2) 

The perturbed metrics now give rise to the second-order perturbations of as 

7 ^ 1 . = -|- ( 57 ^ 1 . -|- -I- ©(third-order perturbations), (5.3) 


where 7 ^ 1 / is the background quantity dehned in eq. (3.4) and satisfies eq. (4.6), and 
hence 5^°‘h = 0 for a,b = 2, 3. The interaction term in the equations of motion for g^i, can 
be calculated explicitly even at second order, and is given by 


X, _ 

^( 9 ) ^ ~ 


/o 0 0 0 \ 

0 0 0 0 


(5.4) 


This tensor satisfies the conditions assumed in appendix B, which, together with the Bianchi 
identity, yield 


= 0 . 

Even at second order, is proportional to 


(5.5) 


7 /)= - 


K?A‘^\ad — hc\ 


X, 

^(9) 


(5.6) 


leading to = 0 as well. These conditions provide the relation between g^‘^'> ah and 

as follows: 


^(2)a^ = - 


- {l^C - A)5^B)Sj^b, 




B 


— A{a + d) + ad — be 


(5.7) 


for a,b = 2,3 and A,B,C = 0,1. Thus, the metric perturbations obey the perturbed Einstein 
equations also at second order, and the number of graviton degrees of freedom coincides with 
that of two massless gravitons even at second order. 

This fact indicates that one cannot discriminate this class of solutions from the cor¬ 
responding solutions in GR even at second order, unfortunately. On the other hand, this 
fact, fortunately, implies that our solutions are free from non-linear instabilities even in cubic 
action, which plague many cosmological solutions in massive gravity, such as the diagonal 
open ELRW solution [16-18], flat FLRW solution [24, 28], and de Sitter solution [23, 28]. 


6 Conclusions and discussion 

In the present study we have investigated the perturbations of a class of spherically symmetric 
solutions in massive gravity and bi-gravity. First, we classified spherically symmetric solutions 
in massive gravity and bi-gravity and identified the specific class for which the background 
equations of motion are identical to a set of the Einstein equations with a cosmological 
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constant. These solutions are allowed only with the one-parameter family of as and 04 
satisfying eqs. (3.11). This class of solutions includes many known solutions, e.g., the FLRW 
solutions in ref. [26, 27], the Schwarzschild(-de Sitter) solutions in ref. [76-79], the LTB 
solution in ref. [27], and the RN solution in ref. [77]. In fact, any spherically symmetric 
solution in GR is included in this class with a suitable choice of the fiducial metric 

Next, we have investigated linear perturbations on this class of solutions. We have 
found that the interaction terms in the equations of motion for both metrics, and 

vanish thanks to the Bianchi identities, and hence the equations of motion reduce 
to eqs. (4.9)-(4.11), which are the perturbed Einstein equations with the relation (4.11). 

We have also found that, in addition to the usual gauge symmetry associated with 
spacetime coordinate transformation, there is another gauge symmetry of the linear pertur¬ 
bations given by eqs. (C. 8 )-(C.ll), which has already been known for the perturbations of 
the Schwarzschild de Sitter solution in dRGT massive gravity [79]. 

We have shown that the above result applies to second-order perturbations as well. 
Thus, one cannot distinguish this class of solutions in massive gravity and bi-gravity from the 
corresponding solutions of GR up to second order. This fact, however, implies that this class 
of solutions do not suffer from the non-linear instabilities, which often appear in the other 
cosmological solutions in massive gravity and bi-gravity. These aspects would suggest that 
massive gravity or bi-gravity with this one parameter family in ( 03 , 04 ) may have additional 
fully non-linear symmetry, which may be responsible for the stability. Further investigations 
are necessary to clarify this point. 

In this article, only spherically symmetric background solutions are discussed. So, it is 
an interesting and open question whether the results obtained in this article hold for more 
general background solutions. Our analysis on the background solutions can at least be ap¬ 
plied to any 7 having the form of eq. (3.4) in any basis vectors, because the background 
equations of motion (3.12) can be obtained in an algebraic way from eq. (3.4) irrespective of 
a concrete expression for 7 . Extending the above analysis to linear and non-linear perturba¬ 
tions of more general solutions, however, is a non-trivial issue simply because such analysis 
accompanies the derivatives. We will address these issues in a future publication. 
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A Concrete examples of background solutions 


A.l Bi-cosmological solutions 

First we consider a family of bi-FLRW solutions, in which physical metric takes the following 
FLRW form: 


= —dt^ + fl^(^) 


dr 


2 


+ r^dn‘^ 


g^udx^dx'^ 


1 - Ar2 


(A.l) 





Comparing this metric with eq. (3.1) yields R{t,r) = a{t)r. We assume that f^i, takes the 
same FLRW metric but in a coordinate {t{t, r), r{t, r), 9, (j)) different from that of 


ff.udx^dx’' 


—dP + P{i) 


dP 

1 - 




foodP + 2foidtdr + fudP + P{i)P{t,r)dQ'^, 


(A.2) 


where 


didi 

Joi = + 


di\^ P{i{t,r)) 


1 — r) 

b‘^{i{t,r)) didi 


dr 


dt dr I — Kf‘^{t,r) dt dr' 

2 , 2 / 


fn — - 


m 

dr 


+ 


P{t{t,r)) i dr\ 


\dr J 


(A.3) 

(A.4) 

(A.5) 


1 — r) 

In order to apply the results of the main body, the radial coordinate r is determined to satisfy 
the following relation, 


AR{t,r) ^ a{t)r 
b{i{t,r)) b{i{t,r))' 


(A.6) 


while the time coordinate t is arbitrary. In this case, the equations of motion for both 
metrics become Einstein equations with cosmological constants so that a{t) and b{t) obey 
the Friedmann equation with respect to each (cosmic) time, t or t. This kind of bi-FLRW 
solution becomes a slight generalization of that found in ref. [36], in which a specific choice 
of the coordinate t is adopted. 

For 6 = 1, A' = 0, and A^^ = 0, the fiducial metric becomes the flat Minkowski 
one and hence this bi-cosmological solution includes that obtained in ref. [26, 27] in dRGT 
massive gravity with the flat fiducial metric. 


A.2 Bi-Schwarzschild de Sitter solutions 

Our results are applied to the following bi-Schwarzschild de Sitter metrics as well: 


g^^dx^dx^ = — 

(l +rW, 

\ r ’« ) i_rM + A(9V2 

(A.7) 

f^iydx^dx'' = - 


(A.8) 


II 

(A.9) 


where and represent Schwarzschild radii, and A^g and A^g are effective cosmological 
constants defined in eqs. (3.13) and (3.15), respectively. Since the Schwarzschild-de Sitter 
metric is a solution of Einstein equation with cosmological constant, this is a vacuum solution 
in our setting. As is the case with the cosmological solution, this black hole solution can 
be obtained with arbitrary choice of the time coordinate t{t,r). By tuning the parameters 
Po and Pi, we can set A^g and A^g to be zeros simultaneously, which corresponds to a 
bi-Schwarzschild solution. 
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B Bianchi identity 


In this appendix, we will show that a symmetric tensor satisfying a condition given below 
must vanish as long as it obeys Bianchi identity and the background metric takes the 
matrix form of eq. (3.1). 

Let us consider the following symmetric tensor 

+ 0(6"+^) (B.l) 

with 

^(n)O^ = = 0, (B.2) 

where e denotes the order of perturbations and A is a constant. The goal of this section is 
to show that vanishes if the Bianchi identity, V^X^,y = 0, is imposed. 

The tensor is symmetric because 

Xp, = gppXP, (B.3) 

= Kgp, + gppX^'^^P, + 0{e^+^), (B.4) 


and both of Xp^, and ^.g^^ are symmetric. Then, from the property of the background metric 
(jpu, it is characterized by three arbitrary functions as follows: 


or equivalently. 




/o 0 0 ^ \ 

000 0 

OOx£^(t,r,0,<^) Xif{t,r,9,(j)) 
\0 0 X^f{t,r,e,4>) Xif{t,r,9,<i>)j 


X(n)/.^ 


/O 0 0 0 \ 

0 0 0 0 


^ ^ R-‘it,r) WW) 

. . X^\t,rfi,<t>) xi'^\t,r,e,4,) 

\ WityiXrXe B?(t,r)sin^e/ 


On the other hand, from the eq. (B.l), the Bianchi identity reads 

V^X'^^ = e^V^xW^^ + 0(e"+^) = 0, 


(B.5) 


(B.6) 


(B.7) 


where is the covariant derivative with respect to g^iy. Then, zero-th and first components 
of this equation are 


V^X^o = -(Xg) + (sin0)-2xg))^ 
V^X^ = -(Xg) + (sin0)-2xg))^ 
which yields the following solution when R is not a constant, 

X^\t,r,9A) 


e" + 0(e”+g = 0, 

(B.8) 

:e" + C)(e”+i) -0, 

(B.9) 


Xg)(f,r,0,0) = -: 


sin"^ 9 


(B.IO) 
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The remaining components of this equation are given by 


= {Rsm6)-^ e" + 0(e”+i) = 0, (B.ll) 

Vf,Xf^3 = {Rsine)-^ + sin 0 ^ 0 (sin 0 X^ 3 ^)) e” + 0(e’"+^) = 0, (B.12) 


where we have used the relation (B.IO). Removing X^'^ from these equations leads to the 

(n) 

following equation for ^33 h 

-^deisinedgX^'^'^) + = 0. (B.13) 

smtl smtl^ 

Since this is just the Laplace equation on a sphere, its solution is constant over the sphere: 

Xg)=/(t,r). (B.14) 

By plugging this solution into eqs. (B.ll) and (B.12), we obtain 


yin) _ 9{t,r) 

^23 - • 

Thus, the solution of the Bianchi identity is given by 


xW = 


/o 0 0 0 \ 

0 0 0 0 

9(t,r) 

Sin*y Sin6» 

yno 


0 0 -44^ 


(B.15) 


(B.16) 


However, the components with sin0 in their denominators are singular at 0 = 0 , 7 r unless 


fit,r) = 0 , 
g{t,r) = 0 . 


Therefore, the regular solution of X^X^y = 0 is 

= 0 


(B.17) 

(B.18) 


(B.19) 


C Additional gauge symmetry of linear perturbations 

The linear perturbations have an additional gauge symmetry, which is combination of gauge 
transformation of g^i, and f^i, separately but keeping the equation (4.11). In this appendix, 
we will give a concrete form of such coordinate transformation. 

For this purpose, let us consider infinitesimal gauge transformation generated by 
for g^j^i, and —)• — (^^ + 6^^) for 

We denote the difference A^5gab — Sfab in eq. (4.11) under this transformation by 
that is. 


A?5gab — Sfab A?5gab — ^fab + Xab- (C.l) 

^To determine the gauge transformation, one establish a bi-tangent bundle T^M i.e. a hbre bundle locally 
isomorphic to M x Tp x Tp'^k To have a usual tangent bundle TM, two horizontal lifts 7r^)(M) and 7r^’^)(M) 
are identified by this relation of the diffeomorphisms so that it determines a diffeomorphism group of the base 
manifold. 
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The additional gauge symmetry is characterized by = 0. The (2, 2) component of this 
condition is given by 


0 = 


i22 


= SfdtR + + Rde6f. 


-2A^R 

The remaining (2,3) and (3, 3) components are given by 


0 = 

0 = 


i23 


—A^R? sin0 
^33 

-2A^R^ siAe 


= d, 


se 


sin 

= -de 


+ sin 6de5^^ 


se \ 

sin 9 ) sin 9 


(C.2) 

(C.3) 

(C.4) 


where we have used eq. (C.2). One can easily find, similarly to eq. (B.13), that these equations 
reduce to the Laplace equation on a sphere: 


1 


sin0 


(sin = 0, 
^ ^ sm^9 


(C.5) 


whose solution becomes 

5e = P{t,r). (C.6) 

Plugging this solution into eqs. (C.3) and (C.4) we find 

6^"^ = Q{t,r)sm9. (C.7) 

To sum up, this additional gauge symmetry is characterized by H(t, r, 9, cj)), P{t, r), Q{t, r) as 


= E{t,r,9,(l)), 


_ dtR{t,r)E{t , r, 9, cj)) + R{t, r)Q{t, r) cos 9 


drR{t, r) 


6^^ = Q{t, r) sin 9, 
Se =P{t,r). 


(C.8) 

(C.9) 

(C.IO) 

(C.ll) 


One may regard R{t,r) itself as a radial coordinate and, in the new coordinates {t,R,9,4>), 
the above transformation (C.8)-(C.ll) with P{t,r) = Q{t,r) = 0 simply reduces to the 
transformation of the time coordinate. 

We can directly observe this symmetry in the action. Actually, the quadratic action of 
the mass term for the linear perturbations becomes 

I ( 0 , 12 ) 

where C{t,r) is the function defined in eq. (3.9), and A^g are effective cosmological 

constants defined in eqs. (3.13) and (3.15), -y/—and a-i'e quadratic perturbations 

of A—g and A~f- Clearly, this term is invariant under the transformation (C.8)-(C.ll) 
because this transformation leaves 6'y°‘b unchanged. 
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D Hamiltonian analysis of linear perturbations 


We will count the number of graviton degrees of freedom of linear perturbations by means of 
the Hamiltonian analysis. So, we omit the matter action in this appendix. For this purpose, 
it is useful to decompose the perturbations in terms of spherical harmonics Yj^ as done in 
ref. [81]. Due to the spherical symmetry of the background metrics, the modes with different 
eigenvalues of rotation {I, m) or parity (odd or even) develop independently, and the dynamics 
of each mode does not depend on m. Hence, we may suppose that m is equal to zero, without 
loss of generality. 


D.l Odd mode perturbations 

Non-vanishing components of the odd mode perturbations with m = 0 are given by 

A9),1 


^903 = ’ {t,r)sm9d0Pi{cosO), 


i>i 


^913 = sm9dePi{cos9), 


i>i 


x ^ • 2 an (dePi{cos9) 

^923 = V h)^ (t, r) sm 9d0 -—- 

\ sm 9 


l>2 


and 


Sfo3 = '^hl/^’^{t,r)sm9d0Pi{cos9), 


i>i 


Sfi 3 = sin 9d0Pi{cos 9), 


i>i 


<5/23 = ^ r) sin^ 9 d 0 


l>2 


d0Pi{cos9) 


sm ( 


(D.l) 

(D.2) 

(D.3) 

(D.4) 

(D.5) 

(D.6) 


where Pi is the Legendre polynomial. In this subsection, hereafter, we omit the suffix I and 
the summation with respect to I for brevity. From the perturbed Einstein-Hilbert action 
with the mass term (C.12), the conjugate momenta of {I = 0,1, 2) are calculated as 


pM-JP-O 

pb) _ _ 2M| 


(2{ln , 


pid) _ 

^2 — 


«« - (g«<>h<<» + - g»'h!f') . 


Sk 


(g) 




py> = 4i=o. 


Sk 


(0 


= ^ = UinpyhiP - h[P' + h^P) 

SPY V-,/ ^ ^ 




_ 6S _ 2XK^M‘f 


p\J ) _ 


pi 


A^P? 


_/ f/oo/,(/) + _ joijyp, 


(D.7) 

(D.8) 

(D.9) 

(D.IO) 

(D.ll) 

(D.12) 
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where A := (Z — 1)(Z + 2), Mpj ;= y/—g represents the determinant of only 

0,1 components: 

:= V-900911 + (D.13) 

We schematically decompose the Hamiltonian density as follows: 

^ “ ^GR,{g) + ^GR,{f) + ^mass) 

where represents the contribution from each Einstein-Hilbert term and the effec¬ 

tive cosmological term, which is the second term (or third term) in the right hand side of 
eq. (C.12). represents the contribution from the hrst term in the right hand side of 

eq. (C.12). This decomposition is justihed because Smass does not include time derivative of 
and /^//. From the expression of the action (C.12), explicitly calculated as 


^mai = MplAC(t,r) 






^ - K 


if) 


(D.15) 


It should be noted that, for Z = 1 mode, vanishes, which implies that dynamics of Z = 1 

mode coincides with that of GR. Therefore, there should be additional gauge symmetry, under 
which each metric transforms independently. This transformation, actually, corresponds to 
the arbitrary function P{t,r) in eq. (C.ll). 

From now on, we focus on Z > 2 modes. The Hamiltonian density from the Einstein- 
Hilbert term is calculated as 


n 


odd 

GR,(g) 




4M2 


v^iprf + 


, (pto))2 + |0i(_/,(9) + 

4M2A 511 1^212 


,(9)/^ dG) 


pi 


511 




511^2' 

\A9)fAg)\2 




R2 


+Ml” + C'c;;; [pj" ,pi’\h'f\h'f], 

where M 2 are some functions of f,r and is given by 


( 9 )^( 1 ) rols) dG) 1 .( 9 ) 1 ( 9)1 


1 "-2 
(D.16) 


"( 9 ) 


CpJ = 2(lnR)'Pf ^ + P^^>' - P^^> - 4M, 


5 ( 9 ) 


( 9 )/ 


yia) 


'ig) 


pi' 


+4AMp2 V^5°^5yl(lnR)^Zi(") - 2Mp2 ^9, h^\ (D.17) 


R2 

with n4 = 0,1. is obtained by replacing g 

The primary constraints of this system are 

/^(o)_ pig) 

^(9) ■“ ^0 

/-a(0)_ pif) 

^if) — -^0 


r2 {R^R) j^ig), 

pV"'- \^J 

>f,R^ AR, Mp2 ^ 

0 , 

’0, 


and then, the total Hamiltonian is 

^odd _ ^odd _|_ 


dr 


(t, r)cf°J -7 (t, r)c{Jj 


^odd _ ^odd I ^odd , rro 

^ — ^GR,(g) ^GR,{f) -^mass; 


'{ 9 ) 

rodd 


U 

rL, 


odd 

GRfg/f) 


= / drR 


odd 

GRfg/f)^ 


ROdd ^ / J RC 

mass / ^ '‘'mass* 


/odd 


(D. 18 ) 

(D.19) 

(D.20) 

(D.21) 

(D.22) 
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Time evolution of the primary constraints is given by 


C 


m 


— rroddi ~ rp'- 


hi) 


( 9 //) 


( 9 //) pig/f) uig/f) l(9//)] 
) ^2 > "'1 ’ "'2 J 


which generate the following two secondary constraints, 


C 


( 1 ) 

( 9 //) 


0 . 


Time evolution of is given by 


^( 1 ) _ 

^( 9 ) “ 


dC 


( 1 ) 


_L /pd) poddi ^ J p(9) rjodd i 


= 2AM, 


' '^'"( 5 ) 
2 \/^ 




and that of c|j.| is given by 


^( 1 ) 

'(f) 


v^, 




(/)^ 


These equations impose another constraint, 

:= ^ 24 ®) - 4-^) « 0 . 

From time evolution of C^‘^\ 

q(2) ^ {c(2)^^odd| 

~ a 2S h(3) podd 1 S h(f) podd i 

/i \n2 ;-«Gi?,(g)i 1^2 ^^GR,{f)J 

4^^ - ^4""^ + 

° 511 ^ 2\Mlgn J 

_ ( *(/) _ ill,4/) + pA 

° hi ^ 2XKm^Jn ^ 


we obtain yet another constraint. 


C(3) := 4^) - + 


VgR" 


-p. 


( 9 ) 


y “ 511 ^ 2AM2 j5ii 

_ ( 4/) _ iiA</) + j/I££ 

“ 111 ‘ 2Ak2MJ/ii 




if) 


(D.23) 

(D.24) 


(D.25) 


(D.26) 


(D.27) 


(D.28) 


(D.29) 


Since includes /iq®^ and terms, the Poisson brackets of and primary constraints 


C 


0 "-0 

(g/f) vanish. Thus, the consistency relation on 


c(3) + {C(3), H} + A^v^g) - v^f) « 0, (D.30) 

determines the combination of multipliers, — '^(f)- Then, no further constraints are 

generated. 
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Since one multiplier remains undetermined, one can easily find that there is gauge 
symmetry in this system. More explicitly, one can confirm that there are two first class 
constraints (and four second class constraints) in this system through the presence of two 
zero eigenvalues of 6 x 6 matrix {C/,Cj}, where Cj represent all of the six constraints. These 
two gauge symmetries correspond to the ones which the theory originally possesses. To 
summarize, the number of graviton degrees of freedom in this system is 


1 

2 



constraints 




(D.31) 


and completely coincides with the case of two massless gravitons. 

For the I = 1 mode, there are four variables (eight variables in phase space), ■f), 

As mentioned above, the interaction term Smass vanishes for I = 1 mode, and hence the action 
reduces to decoupled two Einstein-Hilbert action. Then, there are four first class constraints 
and four gauge symmetries which correspond to the general covariance of g^i, and sep¬ 
arately. These four gauge symmetries can be arranged into the ones of the full theory and 
the additional ones described by P{t,r) in eq.(C.ll). To summarize, the number of degrees 
of freedom of the odd I = 1 mode is 


1 

2 



= 0 . 


(D.32) 


D.2 Even mode perturbations 

Similarly we consider even mode perturbations. Non-vanishing components of the even mode 
perturbations are given by 


6goo = (D.33) 

/>o 

<5501 = ^H[^^\t,r)Piicose), (D.34) 

l>0 

^902 = {t,r)dgPi{cose), (D.35) 

/>i 

^911 = Pi {cos 9), (D.36) 

l>0 

^912 = '^H^^^’\t,r)dBPi{cos9), (D.37) 

i>i 

5922 = Y. r)Pi{cos 9)+ Y dedePi{cos 9), (D.38) 

/>0 l>2 

5933 = 'Y, r) sin^ 9Pi{cos 9) + 'Y, sm9 cos 9 dg Pi {cos 9), (D.39) 

l>0 l>2 


and similar expansions are applied for 5f^u- In this subsection, hereafter, we omit the suffix 
I and the summation with respect to I for brevity. First we treat the I > 2 modes, and the 
Hamiltonian density for even modes is decomposed into 

- ^GR,ig) + ^GR,if) + ttinass- (,U.4Uj 
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'^G^{g/f) I'sprssents a contribution from the Einstein-Hilbert term and effective cosmological 
constant terms, explicitly given by 


^even _ Trid/Dr’i^) 

^GR,{g/f) - t, 


'0 


0 ,( 5 //) 


( 9 //) jjid/f) uio/f) jjis/f) fjid/f) uig/f) 


^pyy/j) pyy/jj p 


I zj\y/J) fry 
> -^5 > ^6 


-H 


( 9 //)^(l) 


1 .( 9 //) 


tpiy/f) pid/f) piy/f) priy/f) rriy/f) priy/f) priy/f)] 
[^3 5 M , r -5 , rig , rig , rig j 


-if, 


( 9 //)^(l) 


tpiy/f) p{y/f) pryy/j) pryy/j) pryy/jj py 

2,(9//) p 4 !-"3 >-"4 >-^6 J 

second order terms of 


riy/f) piy/f) 


riy/f) piy/f)] 


+ 


(D.41) 


with 


{diR)P^^^ + (linear terms of Pg^^^, ) , (D.42) 

= -2Rg^\diR)P^^^ + (linear terms oi P^^\ P^^\ H'f\ hI^\ ,(D.43) 


(1) _ 
1 ,( 9 ) 


C, 


( 1 ) _ 
2 .( 9 ) 


= —2pA^ + (linear terms of pA\ Hn^\ hA\ H, 


( 9 ) p(y) piy) p(y) p{y) 


6 


(D.44) 


where are the conjugate momenta of On the other hand, is given by 

.2 V^C{t, r) 


njeven _ _ _ 

^mass 


A^R^ 


- l{l + 1 ){A^HY’ - - HY’) 




(D.45) 


where Mp = MpK/{l + 21) . Then, the following six primary constraints are imposed. 


^( 0 ) 

^i,{y/f) 

for / = 0,1,2. The total Hamiltonian is 


._ piy/f) 


jjGven _ jje 


/ 


+ dr 


nf^)(t,r)P;®)+nf^)(i,r)Py) 


Time evolution of primary constraints is 


C 


( 0 ) 


= dtC 


( 0 ) 


+ ^^r.(9//)’ 


•'h( 9 //) “ ^^^i,iy/f) 
which impose six secondary constraints, 

C 


Ht} « {Pj 


(y/f) 


Tjeven 

^GR,{g/f) 


}=C 


( 1 ) 

h( 9 //)’ 


( 1 ) 

I,(y/f) 


0 . 


Time evolution of these constraints is given by 

4 ( 1 ) 

" 0 ,( 9 ) 


= -2M^1 f^dsP \- hP) - {A‘^hP - hP) 


^Uy) = -4^pi 


P' HP3 

2 '/~^C{t,r) 


2 17(9) u{f)\ + 1 ) //i 2 [7(9) uU)y 


41) = 2Mpl 


AR^ 

2 V^C{t,r) 


AR^ 


--Y^dBR[{A^HY’ - HY’) - - Pg 

l{l + 1 ) \{A^hP - hP) - {A^hP - hP) 


(D.46) 

(D.47) 

(D.48) 

(D.49) 

(D.50) 

(D.51) 

(D.52) 
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with B = 0,1, and similar terms appear in the constraints for Consequently, we obtain 
two additional constraints: 


:= (D.53) 

cP ■.= (D.54) 

The time evolutions of these constraints are given by 

= A^ + - (^A^Rf^dBRH^^^ + 2A‘^Rf^dBRH[^^^ 

+linear terms of ^ , (D.55) 

4^^ = 2A^H^^'^ - 2H^/^ + linear terms of > (0-56) 

which impose further two constraints, 

^( 3 ) _ ^( 2 ) _ 

Cf) := 4^) « 0. (D.58) 

Since the above constraints include time development of these con¬ 
straints only determines two of the multipliers and hence no more constraint appears. 

One can see that four of the multipliers remain undetermined, which implies that this 

system has corresponding gauge symmetry. Concrete calculation shows that this system has 
eight first class constraints (and eight second class constraints) through eight non-zero eigen¬ 
values of 16 X 16 matrix {C/,Cj}, where Cj represent all of the sixteen constraints. These 
eight constraints are composed of six gauge symmetry of full theory and two additional sym¬ 
metry described by H in eqs. (C.8) and (C.9). To summarize, the number of graviton degrees 
of freedom for even modes can be estimated as 


1 

2 



constraints 




(D.59) 


which again coincides with that of two massless gravitons for I > 2 modes. 

The structure of Hamiltonian analysis is similar for / = 0,1 mode. For I = 0 mode, 
initially we have eight variables (sixteen phase space variables), 


Similar analysis shows that there are ten constraints, 4*^4/) 


C, 


( 1 ) 


C 


( 1 ) 


1,(9//)’ ^0,(g//)’ ^1,(5//) 


o '-'1 


( 2 ) 


(3) 

C\ and six gauge degrees of freedom, that is, there are six first class constraints and four 
second class constraints. Four gauge degrees of freedom come from the ones of the full theory 
and two come from the additional ones described by H in eqs. (C.8),(C.9). Then, the number 
of dynamical degrees of freedom is 


1 

2 



constraints 




(D.60) 


For I = 1 mode, we have twelve variables (twenty-four phase space variables), 


H. 


{ 9 If) jjig/f) ijig/f) u(9/f) 


m 




fourteen constraints, C, 


(0) 


C 


( 0 ) 


C. 


( 0 ) 


C, 


( 1 ) 




)’ '-1, 


( 1 ) 


(9//)’ 
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^ 2 ^ 9 //)’ gauge degrees of freedom, that is, there are ten hrst class constraints 

and four second class constraints. Then, the nnmber of dynamical degrees of freedom is 


1 

2 



constraints 




(D.61) 


It should be noted that six gauge symmetries correspond to the one of full theory, two 
gauge symmetries correspond to H in eqs. (C.8),(C.9), and the other two gauge symmetries 
correspond to Q{t,r) in eq. (C.IO). 
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